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Abstract 

Heat shock factors (Hsfs) play a central regulatory role in acquired thermotolerance. To understand the role of the 
major molecular players in wheat adaptation to heat stress, the Hsf family was investigated in Triticum aestivum. 
Bioinformatic and phylogenetic analyses identified 56 TaHsf members, which are classified into A, B, and C classes. 
Many TaHsfs were constitutively expressed. Subclass A6 members were predominantly expressed in the endosperm 
under non-stress conditions. Upon heat stress, the transcript levels of A2 and A6 members became the dominant Hsfs, 
suggesting an important regulatory role during heat stress. Many TaHsf A members as well as B1 , C1 , and C2 members 
were also up-regulated during drought and salt stresses. The heat-induced expression profiles of many heat shock 
protein (Hsp) genes were paralleled by those of A2 and A6 members. Transactivation analysis revealed that in addi- 
tion to TaHsfA members (A2b and A4e), overexpression of TaHsfC2a activated expression of TaHsp promoter-driven 
reporter genes under non-stress conditions, while TaHsfBIb and TaHsfClb did not. Functional heat shock elements 
(HSEs) interacting with TaHsfA2b were identified in four TaHsp promoters. Promoter mutagenesis analysis demon- 
strated that an atypical HSE (GAACATTTTGGAA) in the TaHsp17 promoter is functional for heat-inducible expression 
and transactivation by Hsf proteins. The transactivation of Hsp promoter-driven reporter genes by TaHsfC2a also 
relied on the presence of HSE. An activation motif in the C-terminal domain of TaHsfC2a was identified by amino 
residue substitution analysis. These data demonstrate the role of HsfA and HsfC2 in regulation of Hsp genes in wheat. 

Keywords: Gene expression, gene regulation, heat shock factors, heat shock proteins, heat stress, transcription factors, 
wheat. 



Introduction 

Temperate cereal crops such as wheat often encounter heat 
stress during the reproductive stage in warm-climate wheat 
production regions (Wardlaw and Wrigley, 1994). Heat stress 
has a significant adverse impact on carbon assimilation and 
starch synthesis in these environments, which leads to reduc- 
tion of grain yield as well as grain quality (Wardlaw and 
Wrigley, 1994; Skylas et al, 2002). At the biochemical level, 
high temperatures cause denaturation of many heat-labile 



proteins and an elevated level of harmful reactive oxygen spe- 
cies (ROS) in plant cells (Xue and Mclntyre, 2011; Mittler 
et al, 2012; Grover et al, 2013). The acquisition of thermo- 
tolerance in plants relies on acclimatization to permissive 
high temperatures, during which time a large number of heat 
protection genes such as those encoding heat shock proteins 
(Hsps) and ROS scavengers are induced (Kotak et al, 2007; 
Mittler et al, 2012). In this heat acclimatization process, heat 
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shock factors (Hsfs) play an important role in regulation of 
this heat-induced transcriptional reprogramming (Kotak 
et al, 2007). 

Hsf proteins are transcription factors, which are present in 
all eukaryotes. Hsfs are characterized by the presence of a Hsf 
DNA-binding domain, which is composed of a central helix- 
turn-helix motif and an adjacent bipartite oligomerization 
domain, made up of of hydrophobic heptad repeats (HR-A 
and HR-B) (Scharf et al. , 2012). The oligomerization domain 
is required for the formation of a trimeric Hsf structure for 
high affinity binding to DNA (Peteranderl and Nelson, 1992; 
Peteranderl et al, 1999). Hsfs are known to recognize the 
multiple inverted repeats of the nGAAn sequence, known 
as the heat shock element (HSE), which appears to be pre- 
sent in the promoters of many heat-inducible genes (Santoro 
et al, 1998; Nover et al, 2001; Guo et ah, 2008; Mittal et ah, 
2011). In general, at least three nGAAn repeats are required 
for effective Hsf binding in Drosophila (Xiao et al, 1991) and 
in Arabidopsis (Kumar et al, 2009). 

The Hsf family in the model plant species, Arabidopsis 
and rice, consists of 21 and 25 members, respectively (Scharf 
et al, 2012), divided into three classes: HsfA, HsfB, and 
HsfC. Among Hsf members, HsfA genes have been func- 
tionally characterized extensively in the model plant species, 
Arabidopsis, serving as bona fide transcription factors bind- 
ing to HSEs. Despite the DNA-binding domain of the three 
classes of Hsf members being highly conserved, knowledge 
about the biological role of the HsfB and HsfC members is still 
scarce. It is generally recognized that HsfB and HsfC members 
do not contain an aromatic (W, Y, and F)-hydrophobic (L, I, 
V, and M)-acidic (D and E) residue-rich (AHA) activation 
domain present in HsfA members (Kotak et al, 2004). HsfB 
members have been demonstrated to act as repressors (Ikeda 
and Ohme-Takagi, 2009; Ikeda et al, 2011; Zhu et al, 2012). 
Recently, Mittal et al. (2011) have shown that HsfB4b from 
rice is also capable of binding to HSE. The binding sequence 
of Arabidopsis HsfBl has been shown to be different from 
that of HSEs; it binds to the TL1 element (GAAGAAGAA) 
and regulates pathogen defence genes (Pajerowska-Mukhtar 
et al, 2012; Pick et al, 2012). However, Arabidopsis HsfBl 
can positively regulate the acquired thermotolerance (Ikeda 
et al, 2011). 

Some Hsf genes are constitutively expressed. A number 
of constitutively expressed HsfAl proteins such as Hsf Ala 
in tomato are known to be maintained in an inactive mono- 
mer state by association with Hsp90/Hsp70 under non-heat 
stress conditions (Hahn et al, 2011). This complex is dis- 
sociated upon heat stress, which enables HsfAl a to form an 
active trimer. In Arabidopsis, constitutively expressed Hsf Ala 
and Hsf Alb are known to be the early response Hsf genes 
(Lohmann et al, 2004; Busch et al, 2005), serving as tran- 
scriptional activators for high-level induction of late response 
Hsf members such as HsfAl, which also plays an important 
role in acquired thermotolerance (Charng et al, 2007). 

In plants, there are some species-specific features in the role 
of individual Hsf members in regulating genes involved in the 
heat stress response. HsfAl a in tomato has been described 
as a master regulator for triggering the heat response and 



acquired thermotolerance (Mishra et al, 2002). However, 
no master regulator has been found in Arabidopsis (Scharf 
et al, 2012). Instead, three subclass Al Hsfs (Ala, Alb, and 
Aid) in Arabidopsis are functionally redundant for triggering 
the heat response (Liu et al, 2011; Liu and Charng, 2012). 
In tomato, HsfBl functions as a transcription co-activator, 
cooperating with class A Hsfs during heat stress (Bharti et al, 
2004), while Arabidopsis HsfBl was found to act as a repres- 
sor (Czarnecka-Verner et al , 2004; Ikeda and Ohme-Takagi, 
2009; Ikeda et al, 2011; Zhu et al, 2012). In addition, the 
compositions of Hsf classes differ markedly among plant 
species (Scharf et al, 2012). Most notably, subclass C2 Hsfs 
appear to be monocot specific. However, the role of this sub- 
class of Hsf in monocots is currently unknown. Conversely, 
subclasses A9, B3, and B5 are not present in monocot species 
(Scharf et al, 2012). HsfA9 in sunflower and Arabidopsis is 
expressed specifically in seeds and controls Hsp expression 
during seed development (Almoguera et al, 2002; Kotak 
et al, 2007). In rice, HsfA7 appears to be expressed specifi- 
cally in the seeds under normal conditions (Chauhan et al, 
2011a). 

Investigation into the molecular basis of heat tolerance in 
wheat was pursued at the Hsp level as early as in the 1980s 
(Zivy, 1987; Krishnan et al, 1989). Acquired thermotol- 
erance in some thermotolerant wheat genotypes appears 
to be linked to higher transcript or protein levels of some 
Hsps (Vierling and Nguyen, 1992; Joshi et al, 1997; Skylas 
et al, 2002). More recently, Chauhan et al (2012) have 
shown that overexpression of a small wheat chloroplastic 
Hsp (TaHSP26) in Arabidopsis improves thermotolerance. 
Comparative expression analysis between two genotypes 
with contrasting thermotolerance using the Affymetrix wheat 
genome array has shown that a large number of probe sets are 
differentially expressed in the heat-stressed leaves between the 
thermotolerant and thermosusceptible genotypes (Qin et al, 
2008). However, to date, no study has reported on the role 
of Hsf genes in regulating Hsp expression in wheat. Little is 
known about the Hsf family structure in wheat and the role 
of individual members in heat stress adaptation. Bread wheat 
(Triticum aestivum) is a hexaploid species and the Hsf fam- 
ily structure in wheat is expected to be more complex than 
that in diploid species. Only two members, HsfA4a and Hsf3 
(HsfB2a), of the wheat Hsf family have been functionally 
analysed. TaHsfA4a is a cadmium-up-regulated gene and is 
involved in cadmium tolerance in wheat (Shim et al, 2009). 
Overexpression of TaHsf3 in transgenic Arabidopsis leads to 
enhanced thermotolerance and freezing tolerance (Zhang 
et al, 2013). 

In this study, a genome-wide identification of the expressed 
members of the wheat Hsf family was performed using a bio- 
informatics analysis of currently available wheat sequences, 
including those isolated from the authors' laboratory. A total 
of 56 members of the Hsf family from T. aestivum were iden- 
tified. The wheat Hsf family structure was determined by phy- 
logenetic analysis in accordance with the Arabidopsis and rice 
Hsf families. In order to understand their potential roles, the 
expression profiles of these Hsf members were investigated in 
wheat in various organs, in response to heat stress and other 
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abiotic stresses (drought and salt stress). Functional analyses 
of representative wheat Hsf members were also performed 
to determine the functional HSE sequences in the promoter 
regions of four Hsp genes, to investigate whether they are 
capable of transactivating the expression of Hsp genes, and 
to identify with which elements they interact. Of particular 
note, it was demonstrated that a HsfC member was a tran- 
scriptional activator of Hsp genes; this has not been shown 
in any other plant species. These analyses provide some fun- 
damental insights into the role of Hsfs in mediating wheat 
adaptation to heat stress through regulation of Hsp genes as 
well as their potential roles in adaptation to other major abi- 
otic stresses. 



Materials and methods 

Plant materials and growth conditions 

Spring wheat (T. aestivum L. cv. Bobwhite) plants were grown in a 
controlled-environment growth room in 1.5 litre pots, containing 
a 3:1:1 mix of sand:soil:peat under night/day conditions of 16n 
light (500 umol ra 2 s _1 ), 16/20 °C, and 80/60% relative humidity. 
For heat treatment experiments, Diatomite [silica gel-based arti- 
ficial soil (0.5-2mm), Burnell Agencies, Brisbane, Australia] was 
used for plant growth for facilitating isolation of clean roots. Heat 
treatment of 1 -month-old plants at 36 °C commenced at 2h after 
lights on, and the leaves and roots of heat-treated plants were har- 
vested at 1.5 h, 5h, and 3 d, with 5h of heat treatment per day. 
Control plants at 20 °C were harvested at the same time as for 
1.5h heat-treated samples. All harvested samples were immedi- 
ately immersed in liquid nitrogen and stored at -80 °C prior to 
RNA extraction. 



Identification and classification of Hsf members from T. aestivum 

Fhe sequences of Hsf DNA-binding domains from rice and 
Arabidopsis Hsf proteins in the Plant Transcription Factor Database 
(Riano-Pachon et al. , 2007) were used for searching T. aestivum Hsf 
expressed sequence tags (ESTs) from the NBCI ESF database and 
the assembled tentative consensus (TC) sequences from the Triticum 
aestivum Gene Index (FaGI) database (http://compbio.dfci. harvard, 
edu/tgi/). Fo increase the coverage of partial wheat Hsf cDNAs that 
lack sequences in the Hsf DNA-binding domain, full-length bar- 
ley Hsf protein sequences retrieved from the NCBI sequence data- 
base were used to find additional wheat Hsf ESFs in the NCBI ESF 
database. Unique Hsf representatives in wheat were identified by 
sequence alignment and contig assembly. After elimination of unre- 
liable sequences at the ends of some poor quality ESFs, sequences 
that had >99.5% nucleotide identities were used for assembly into 
contigs. Fhis analysis corrected many wrongly assembled Hsf FC 
sequences in the FaGI database. Fhe assembled sequences were 
checked for correctness in the wheat genomic sequence database 
(Wilkinson et al, 2012), and partial cDNAs were further extended 
with wheat genomic sequences. Many ESF singletons which had 
no match with the wheat genomic sequence in the CerealsDB 
were discarded, as these singletons most probably represent poor 
sequence quality ESFs. Fo improve the coverage of the full-length 
Hsf cDNAs and the accuracy of Hsf cDNA assembly, the coding 
region sequences of 15 Hsf genes were also isolated using reverse 
transcription-PCR (RF-PCR). Fhese cDNA sequences were then 
incorporated into the Hsf gene assembly. At the final stage, the 
wheat Hsf protein sequences from the assembled Hsf nucleotide 
sequences were checked by BLASFP research in the NCBI protein 
database for the presence of the Hsf DNA-binding domain and the 
partial Hsf sequences were checked for close homologues of barley 
full-length Hsf proteins. 



Phylogenetic analysis was used to classify FaHsf members into 
three classes and further subclasses, based on the rice and Arabidopsis 
Hsf protein classification used by Scharf et al. (2012). Hsf DNA- 
binding domain and heptad repeat region (HR-A/B) sequences of 
Hsf proteins were used for generation of a phylogenetic tree by 
ClustalW alignment and the unrooted Neighbor-Joining method 
using MEGA 5.10 (Famura et al, 2011). Neighbor- Joining analy- 
sis was performed with pairwise deletion and Poisson correction. 
Bootstrap analysis was performed with 1000 replicates to assess the 
level of statistical support for each tree node. Partial FaHsf proteins 
that lack the DNA-binding domain and HR-A/B region were classi- 
fied based on their highest sequence homology to full-length FaHsf 
members, but were excluded from the phylogenetic analysis. 



Total RNA extraction 

Frozen samples were ground to fine powder in liquid nitrogen. 
Fotal RNA was isolated from samples using Plant RNA Reagent 
(Invitrogen, Carlsbad, CA, USA), according to the manufacturer's 
instructions. RNA was further purified through a Qiagen RNeasy 
column (Qiagen, Australia) after pre-treatment with RNase-free 
DNase I (Xue and Loveridge, 2004). 

Isolation of FaHsf cDNAs 

Wheat cDNAs were synthesized from total RNA prepared from the 
leaves and roots of wheat with or without heat treatment. Wheat 
Hsf cDNAs were isolated using RF-PCR or 3 '-RACE (rapid ampli- 
fication of cDNA ends) followed by PCR using primers designed 
from assembled Hsf sequences. Fhe PCR-amplified products were 
cloned into pGEM-F Easy vector (Promega) and sequenced using a 
BigDye terminator cycle sequencing kit (Applied Biosystems, Foster 
City, CA, USA). 

Expression analysis using quantitative real-time PCR 

Fhe transcript levels of wheat genes were quantified from cDNA 
samples synthesized from DNase I-treated total RNA using real- 
time PCR with a ViiA™ 7 system (Applied Biosystems) and SYBR 
Green PCR Master Mix (Applied Biosystems) according to the 
manufacturer's instructions. Fhe sequences of primer pairs used for 
real-time PCR are listed in Supplementary Fable SI available at JXB 
online. Fhe gene-specific primers were designed at the C-terminal 
domain or 3 '-untranslated region. Fhe gene specificity of primers 
for each gene during primer design was checked by blasting primer 
sequences in the FaGI database (http://compbio.dfci.harvard.edu/ 
cgi-bin/tgi/gimain.pl?gudb=wheat) using an expect value setting at 
10 000 and was also checked with all FaHsf DNA sequences assem- 
bled from this study to avoid matching to non-targeted genes. 

Fwo wheat genes (TaRPII36, RNA polymerase II 36kDa subu- 
nit; and TaRP15, RNA polymerase I, II, and III, 15kDa subunit) 
were selected as internal reference genes for calculation of relative 
transcript levels of the genes under study (Xue et al, 2006, 2008a). 
Fhe mRNA levels of these internal reference genes were similar in 
the control and heat-treated samples, as checked by the use of an 
external reference mRNA (626 nucleotides) in vitro transcribed from 
a bovine cDNA (CF767388), which was added to each RNA sample 
before cDNA synthesis (Xue et al, 2011). Fhe PCR efficiency of 
each primer pair was determined by a dilution series of samples. 
Fhe specificity of real-time PCR amplification was confirmed by a 
single peak in melting temperature curve analysis of real-time PCR- 
amplified products. Fhe apparent expression level of each gene rela- 
tive to an internal reference gene (TaRP15) was calculated using the 
following formula: Er ct /Et ct xF (Stephenson et al, 2007), where 
Ct is the cycle threshold (PCR cycle number of reference or target 
gene required for reaching the signal point used for detection across 
samples), Er is the reference gene (TaRPIS) amplification efficiency, 
Et is the target gene amplification efficiency, and F is an amplicon 
size factor (reference gene amplicon size/target gene amplicon size). 
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Apparent expression level values were tentatively used here to pro- 
vide an approximate estimation of relative expression levels among 
various genes under the situation where the absolute quantification 
of mRNA levels for a large number of genes using a cRNA (or 
cDNA) calibration curve is not possible. 

Plasmid construction 

pTaHsfA2b-CELD was constructed by translational fusion of the 
coding region sequence of TaHsf A2b to the N-terminus of the 
6xHis-tagged CELD reporter under the control of the tac pro- 
moter (Xue, 2005). CelD encodes a l,4-|3-glucanase (cellulase) 
from Neocallimastix patriciarum (Xue et al, 1992). pUbiTaHs- 
fA2b, pUbiTaHsfA4e, pUbiTaHsfBlb, pUbiTaHsfClb, and pUbi- 
TaHsfC2a plasmids were constructed by replacing xylanase in 
pUbiSXR (Vickers et al., 2003) with the coding region of TaHsf 
cDNAs. TaHsfC2a mutant constructs (pUbiTaHsfC2a-mCl and 
pUbiTaHsfCa2-mC2) were constructed using PCR and oligonucle- 
otide primers containing the desired sequence mutation. TaHspl 7 
and TaHsp90.1-Al promoter-driven gfp (green fluorescent protein) 
reporters were constructed by replacing the barley HVAl promoter 
in the HVAl promoter-driven gfp reporter gene (Xiao and Xue, 2001) 
with the PCR-amplified fragment of TaHspl 7 or TaHsp90.1-Al 
promoters, isolated in this study (see the Hsp promoter isolation 
section below). TaHsp90.1-A\ promoter mutants (psHsp90gfp and 
pAHsp90gfp) were constructed using PCR-directed promoter trun- 
cation. The pHSE90gfp and pHSE17gfp constructs were made by 
adding three repeats of TaHsp90.1El orTaHspl7El elements to the 
TaHsp90.1-Al minimal promoter construct (pAHsp90gfp). 



DNA-binding activity assays 

For DNA-binding assays 6xHis-tagged TaHsfA2b-CELD fusion 
protein was purified using Ni-NTA magnetic agarose beads 
(Qiagen, Australia) as described previously (Xue, 2005). Biotin- 
labelled double-stranded oligonucleotide probes were synthesized 
by filling in partially double-stranded oligonucleotides using the 
Tag polymerase reaction as follows: a 50 ul PCR containing lx 
Tag polymerase buffer, 200 uM dNTPs, 150 pmol single-stranded 
oligonucleotide [5 '-(sequence containing HSE or its mutants)- 
GAGTGGTGATTCCGGGCCTT (20 bp linker sequence)], 150 
pmol antisense biotin-labelled primer (5'-biotin-AAGGCCCG 
GAATCACCACTC) that was reverse complementary to the linker 
sequence, and 1 U of Tag polymerase (Roche) with PCR cycling 
parameters as described previously (Xue et al., 2006). 

The DNA-binding activity of TaHsfA2b-CELD was measured 
as described previously (Xue, 2002a) using StreptaWell High Bind 
(streptavidin-coated 96-well plates from Roche) and binding/wash- 
ing buffer [25 mM HEPES/KOH, pH 7.0, 50 mM KC1, 2mM MgCl 2 , 
and 0.5 mM dithiothreitol (DTT) containing 0.15 (xg uT 1 shared 
herring sperm DNA, 0.3 mg ml -1 bovine serum albumin, 10% glyc- 
erol, and 0.025% NP-40. About 50000 fluorescent units ir 1 of the 
CELD activity of TaHsf A2b-CELD fusion protein and 0.4 pmol of 
biotinylated probes were used per assay. The cellulase activity of the 
CELD fusion proteins bound to immobilized biotinylated probes 
was assayed by incubation in 100 ul of the CELD substrate solution 
(ImM methylumbelliferyl-P-D-cellobioside in 50 mM Na-citrate 
buffer, pH 6.0) at 40 °C for 4h. DNA-binding assays with a biotin- 
labelled double-stranded oligonucleotide without a HSE were used 
as a control of background activity. 

Transactivation assays 

Transactivation assays were performed as described previously 
(Xue, 2003). Constructs were introduced into the seedlings of wheat 
(cv. Bobwhite) by particle bombardment (Xue et al, 2003). The 
effector gene was co-introduced with a gfp reporter gene driven by 
TaHspl 7, TaHsp90.1-Al, or its derivative promoter to determine 
the transactivation activity. The gfp reporter gene without a TaHsf 



effector gene was used as a control. A |3-glucuronidase (GUS+) 
reporter driven by the maize Ubil promoter, as described by Vickers 
et al. (2003), was also co-bombarded for validation of transforma- 
tion events among assays. The bombarded seedlings were kept at 
room temperature (22 °C) or 36 °C in the dark for ~20h. GFP foci 
were examined under a fluorescence microscope (Xue, 20026). The 
tissue sections that had GFP foci were subsequently stained for his- 
tochemical detection of GUS activity according to the method of 
Jefferson (1987). 

Isolation of TaHSP promoter sequences 

The promoter sequences of TaHspl7 and TaHsp90.1-Al were iso- 
lated with PCR amplification of genomic DNA of T. aestivum geno- 
type SB169 (Xue et al, 20086). PCR primers were designed based 
on the assembled sequences through extension of the TaHspl7 and 
TaHsp90.1-Al EST or cDNA sequences using the wheat genome 
sequence database in CerealDB (Wilkinson et al, 2012). The PCR- 
amplified DNA fragments were cloned and sequenced. Two pro- 
moter sequences were deposited in GenBank [TaHspl 7 promoter 
(1214bp upstream of the translation start codon), KF208539; and 
TaHsp90.1-Al promoter (1375 bp upstream of the translation start 
codon), KF208540]. 

Expression analysis of Affymetrix GeneChip wheat genome 
array data sets 

The wheat genome array (Affymetrix GeneChip) contains 61 
127 probe sets representing 55 052 transcripts for all 42 chromo- 
somes in the wheat genome. The raw GeneChip data of TA23 and 
E-MEXP-971 were retrieved from http://www.plexdb.org and http:// 
www.ebi.ac.uk/arrayexpress/ as eel files, respectively. The Affymetrix 
data sets were analysed as described previously (Xue et al, 2013), 
using GeneChip robust multiarray average (GC-RMA) for nor- 
malization (Wu et al , 2004). Normalized values were converted to 
non-log values for comparative analysis of TaHsf genes between 
treatments. Probe sets with normalized hybridization signals of <20 
were considered as not detectable. 



Results 

Fifty-six members were identified from the TaHsf family 

Fifty-six TaHsf members were identified through a combina- 
tion of sequence analysis of wheat EST and nucleotide col- 
lection databases in NCBI (http://blast.ncbi.nlm.nih.gov/), 
followed by confirmation and sequence extension analyses 
using wheat genome sequences in the CerealDB (http://www. 
cerealsdb.uk.net/CerealsDB/Documents/DOC_search_reads. 
php, Wilkinson et al,, 2012), and isolation of 15 TaHsf 
cDNAs in this study (Table 1). Sequence IDs used for assem- 
bly of these TaHsf members are listed in Supplementary 
Table S2 at JXB online. The deduced protein sequences of 56 
TaHsf genes are shown in Supplementary Fig. SI. Of these 
56 genes, 39 TaHsf genes contain a full-length coding region 
sequence. A number of these partial TaHsf proteins identi- 
fied contain half a Hsf DNA-binding domain. There is a very 
large intron between the LNTY and GFRK sequence in the 
DNA-binding domain, and this makes successful extension 
of partial ESTs difficult using the incomplete wheat genome 
sequence database in the CerealDB. It is possible that two 
partial cDNAs, TaHsfA2i and TaHsfA6f, containing an 
LNTY half-sequence may be a part of GFRK-half cDNAs 
(Supplementary Fig. S2). The multiple sequence alignments 
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Table 1 . The list of TaHsf members identified 



Gene 


Full or partial 


NO. OT 


P' 


IVIOI. 


Affymetrix probe set ID 




tC nr N) DRF 


dllllllU dUIUo 




wt fkna) 


Inn nf nmhp mat/*hoc\ a 

^IIU. UI |JIUUc 1 1 Id IOI IcoJ 


TaHsfAl a 


Full 


529 


4.99 


58.1 




TaHsfAl b 


Full 


522 


4.91 


57.5 


Ta.8059.1.S1_at(7) 


TaHsfAl c 


Almost full 


NA 


NA 


NA 




TaHsfA2a 


Full 


353 


5.36 


39.7 


TaAffx.92707.2.S1_at (11) 


TaHsfA2b 


Full 


413 


4.99 


45.6 


Ta.6737.2.S1_a_at (1 1), TaAffx.105519.1 .S1_at (8) 


TaHsfA2c 


Full 


404 


5.06 


44.8 


TaAffx.105519.1.S1_at(8) 


TaHsfA2d 


Partial (C) 


NA 


NA 


NA 


TaAffx.92707.1.S1_at(9) 


TaHsfA2e 


Partial (C) 


NA 


NA 


NA 


Ta.6737.2.S1_a_at (11) 


TaHsfA2f 


Partial (C) 


NA 


NA 


NA 


Ta.1276.1.A1_at(10) 


TaHsfA2g 


Partial (N) 


NA 


NA 


NA 




TaHsfA2h 


Full 


405 


5.06 


44.9 


TaAffx.105519.1.S1_at(8) 


TaHsfA2i 


Partial (N) 


NA 


NA 


NA 




TaHsfA3a 


Full 


499 


5.7 


54.8 


TaAffx.12296.1.S1_at(11) 


TaHsfA3b 


Partial (C) 


NA 


NA 


NA 




TaHsfA4a 


Full 


432 


5.36 


48.4 




TaHsfA4b 


Full 


441 


5.18 


49.5 




TaHsfA4c 


Full 


448 


4.91 


50.2 


TaAffx.1 20360. 1.A1_at (4) 


TaHsfA4d 


Full 


442 


5.11 


49.7 


TaAffx.5995.1 .S1_at (8), TaAffx.1 20360.1 .A1_at (5) 


TaHsfA4e 


Full 


445 


4.93 


49.9 


TaAffx.1 20360. 1.A1_at (11) 


TaHsfA4f 


Partial (C) 


NA 


NA 


NA 




TaHsfA5a 


Full 


458 


5.21 


49.9 




TaHsfA5b 


Full 


455 


5.33 


49.9 




TaHsfA6a 


Full 


341 


5.02 


39.5 


Ta.27873.1.A1_at (9) 


TaHsfA6b 


Partial (C) 


NA 


NA 


NA 




TaHsfA6c 


Partial (C) 


NA 


NA 


NA 


Ta.28772.1.S1_at (6) 


TaHsfA6d 


Partial (C) 


NA 


NA 


NA 


Ta.28772.1.S1_at(9) 


TaHsfA6e 


Full 


370 


5.03 


42 


Ta.28772.1.S1_at (6) 


TaHsfA6f 


Partial (N) 


NA 


NA 


NA 




TaHsfA7a 


Full 


353 


4.79 


38.2 




TaHsfA7b 


Full 


351 


4.83 


38.2 




TaHsfASa 


Full 


384 


5.15 


43.1 


TaAffx.40092.1.S1_at(7) 


TaHsfASb 


Partial (C) 


NA 


NA 


NA 


TaAffx.40092.1.S1_at(5) 


TaHsfA8c 


Partial (C) 


NA 


NA 


NA 


TaAffx.40092.1.S1_at(11) 


TaHsfBIa 


Full 


298 


9.5 


32.2 


Ta.11671.1.S1_at(9) 


TaHsfBI b 


Full 


298 


9.31 


32.3 




TaHsfBIc 


Full 


298 


9.2 


32.1 


Ta.11671.2.S1_x_at(6) 


TaHsfB2a 


Full 


314 


9.46 


34.7 




TaHsfB2b 


Full 


374 


5.44 


40.4 




TaHsfB2c 


Full 


397 


4.89 


41.4 




TaHsfB2d 


Full 


396 


4.89 


41.1 




TaHsfB2e 


Full 


397 


5.28 


41.3 


Ta.18067.2.S1_x_at(4) 


TaHsfB4a 


Full 


320 


6.55 


35.3 




TaHsfB4b 


Partial (N) 


NA 


NA 


NA 




TaHsfB4c 


Partial (N) 


NA 


NA 


NA 




TaHsfCIa 


Full 


325 


6.08 


35.8 




TaHsfClb 


Full 


236 


6.91 


26.1 


Ta.8266.1.A1_s_at (11), TaAffx.34778.1 .S1_at (11) 


TaHsfCIc 


Full 


234 


7.65 


25.8 


Ta.8266.1.A1_s_at(7) 


TaHsfCId 


Full 


241 


8.76 


26.4 


Ta.8266.1.A1_s_at(11) 


TaHsfCle 


Full 


325 


5.94 


35.7 




TaHsfC2a 


Full 


263 


6.11 


28.0 


Ta.5852.1 .A1_at (9), Ta.5852.1 .A1_x_at (9) 


TaHsfC2b 


Full 


270 


6.92 


28.7 




TaHsfC2c 


Almost full 


NA 


NA 


NA 




TaHsfC2d 


Full 


274 


6.98 


30 


Ta.13964.1.S1_at(11) 


TaHsfC2e 


Full 


276 


6.46 


30.3 




TaHsfC2f 


Full 


264 


5.43 


29.1 


Ta.13964.1.S1_at (11) 


TaHsfC2g 


Full 


275 


8.4 


30.2 


Ta.25627.1 .A1_at (1 1), Ta.25627.1 .A1_x_at (11) 



N, the N-terminal part; C, the C-terminal part; NA, not applicable. 

a The Affymetrix wheat genome array has 1 1 probes for each probe set ID. 
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of the DNA-binding domains and heptad repeats (HR-A core 
and HR-B) of TaHsf proteins are shown in Supplementary 
Figs S2 and S3, respectively, illustrating highly conserved 
amino residues in the DNA-binding domain and primary 
sequence features in the HR-A core/B region in these TaHsf 
proteins. 

Phylogenetic analysis and classification of TaHsf 
proteins 

In order to gain some insight into the potential function of 
TaHsf proteins from known function Hsfs in model plant 
species (Arabidopsis and rice), phylogenetic analysis of Hsf 
proteins from wheat, rice, and Arabidopsis was performed 
using the sequence region from the start of DNA-binding 
domain to the end of the HR-A/B region (Fig. 1). Class and 
subclass classification of TaHsf proteins is based on those of 



Arabidopsis and rice Hsf proteins as reported by Scharf et al. 
(2012). The phylogenetic tree of the TaHsf family members is 
shown in Supplementary Fig. S4 at JXB online. From Fig. 1, 
it can be seen that while class B and C Hsf members form dis- 
tinct groups, class A Hsf protein members form at least three 
smaller clusters; two of the class A clusters appear to be more 
similar to class C Hsf proteins than to the third class A clus- 
ter. Class C Hsf proteins appear to be more similar to class 
A Hsfs, compared with class B members. Wheat Hsf proteins 
cluster more closely with rice Hsfs than with Arabidopsis Hsfs. 
Arabidopsis A2 is not clustered with the monocot A2 group. 
Wheat Hsf A, B, and C classes contain eight, three, and two 
subclasses, respectively, which are the same as those in the 
rice Hsf family. Subclasses A9 and B3 are absent in wheat 
and rice, and subclass C2 is unique in monocot species. TaHsf 
members with an incomplete sequence of the DNA-binding 
domain and HR-A/B region were subsequently classified 
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Fig. 1. The Neighbor-Joining phylogenetic tree of Hsf proteins from wheat, rice, and Arabidopsis. The N-proximal regions (from the 
start of the DNA-binding domain to the end of the HR-A/B region) of Hsf proteins were used for construction of the phylogenetic 
tree using the MEGA 5.10 program. Unrooted Neighbor-Joining analysis was performed with pairwise deletion and Poisson 
correction. For rice (prefixed by Os) and Arabidopsis (prefixed by AT) Hsf proteins, both locus ID and subclass number were used (e.g. 
AT4G1 7750A1 a=AtHsfA1 a with a locus ID AT4G1 7750). TaHsf proteins are marked with asterisks. Bootstrap values >50 are shown. 
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based on their closest sequence homology with TaHsf mem- 
bers containing a complete sequence. In comparison with the 
Arabidopsis Hsf family, the wheat Hsf family contains a large 
number of A2 group members (Table 1), in addition to the 
differences in the A9, B3, and C2 groups. 

Organ distribution of TaHsf transcripts 

In order to understand the biological roles of individual 
TaHsf members in wheat, their relative transcript distribu- 
tion was surveyed under non-stress conditions in six major 
organs (leaf, root, stem, hull, endosperm, and embryo), 
with leaf and endosperm studied at two developmental 
stages. Real-time PCR primers targeted to the 3 ' regions or 
C-terminal domain region of cDNAs for TaHsf members 
were designed with full-length cDNA or containing full or 
partial C-terminal part sequences, with the exception of 
TaHsf A2h, which was isolated after completion of expression 
analysis. The mRNA levels of individual members are pre- 
sented as apparent expression levels relative to the internal 
control gene, TaRP15 (Fig. 2), which are used for the assess- 
ment of relative expression levels between genes (Stephenson 
et al, 2007). TaHsf subclass Al group genes were constitu- 
tively expressed at relatively high levels among the six organs. 
Other class A genes constitutively expressed at relatively high 




Young leaf IMo 
Flag leaf 20DAA 
Hull 20DAA 
Young root IMo 
Stem 20DAA 
Endosperm 20DAA 
Endosperm 27DAA 
EmbrYO 42DAA 



Ala Alb Ale A2a A2b A2c A2d A2e A2f A3a A3b A4a A4b A4c 




Bla Bib Blc B2a B2b B2c B2d B2e B4a B4b B4c Cla Clb Clc Cld C2a C2b C2c C2d C2e C2f C2g 

Fig. 2. Relative mRNA abundance of TaHsf members in wheat 
organs. Values are means ±SD of three biological replicates and 
are expressed as apparent expression levels relative to a control 
gene TaRP15. 1 Mo, 1 -month-old; 20DAA, 20 d after anthesis; 
27DAA, 27 d after anthesis; 42DAA, 42 d after anthesis. 



levels in the green organs examined (leaf, stem, and hull) were 
some of the A2, A6, and A8 members. Class A members 
that were predominantly expressed in the endosperm were 
A2b, A2c, A2e, A5b, A6c, A6d, and A6e, and the apparent 
expression levels of many of these endosperm-predominantly 
expressed members were higher than that of the Al members 
in the endosperm. In particular, the constitutive expression 
of A6e in the endosperm was at least three times higher than 
that of the Al group genes. A6e is highly homologous to a 
partial TaHsf cDNA (GD186916) isolated from the T. aesti- 
vum developing seed heat stress forward subtractive library, 
which has been shown to have a seed-predominant expression 
pattern (Chauhan et al, 201 lb). Two A2 members (A2c and 
A2e) were expressed in the maturing embryo at a higher level 
than the Al members. 

Some class B members (B2c, B2d, and B2e) were also con- 
stitutively expressed in these six organs at a relatively high 
level (Fig. 2). The subclass Bl members were expressed 
at higher levels in organs at the reproductive stage than in 
young leaves and roots. The expression of B4 members was 
mainly confined to the roots and embryos under normal con- 
ditions. The transcript distribution of class C members in 
these organs had the following features: (i) several CI mem- 
bers (Clb, Clc, and Cld) were highly expressed in the young 
leaves and embryos, with the apparent expression level of Clc 
in the young leaves higher than those of the Al group; and (ii) 
C2 members were predominantly expressed in embryos, and 
their apparent expression levels in embryos and roots were 
much higher than any other subclasses of TaHsf genes. 

To understand the functionality of TaHsf proteins under 
normal conditions, some representative Hsp genes were 
also analysed. Although many class A members were con- 
stitutively expressed in vegetative organs, many Hsp genes 
(TaHspl6.9, TaHspl7, TaHsp26.6, and TaHsp90.1-Al) were 
not expressed or were expressed at very low levels in these veg- 
etative organs under normal conditions (Fig. 3A). However, 
many of these wheat Hsp genes, particularly TaHsp70d, were 
constitutively expressed at a high level in endosperms and 
maturing embryos. 

Subclasses A2, A6, and B2 members become 
predominant TaHsf transcripts during heat stress 

To provide a heat response profile for individual TaHsf mem- 
bers, the heat responsiveness of these genes in the leaves and 
roots of 1 -month-old plants with short- (1.5 h and 5h) and 
long-term (3 d with 5 h heat treatment each day) heat stress 
(at 36 °C) was examined. The expression of subclass Al group 
members was not up-regulated during heat stress; in fact, the 
Ala mRNA level was down-regulated in leaves with a short- 
term heat treatment and a marked reduction in its mRNA 
level was observed in the 1.5 h treatment (Fig. 4). Transcript 
down-regulation by heat in the leaves was also observed in the 
A3 group, A4c, A4d, A8a, and A8c genes. Heat up-regulation 
was observed in all A2, A6, and A7 members, as well as in 
some members of the A4 and A5 groups. In the early heat 
treatments, the A2 and A6 groups were the predominantly 
expressed HsfAs in the leaves. In the roots of heat-stressed 
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Fig. 3. Expression profiles of TaHsp genes. Values are means 
±SD of three biological replicates and are expressed as apparent 
expression levels relative to a control gene TaRP15. (A) Organ 
imRNA distribution. (B) Response to heat stress in the leaves and 
roots. Statistical significance in differences between control and 
heat-treated groups (36 °C for 1 .5h, 5h, or 3 d with 5h d~ 1 ) is 
indicated by an asterisk. 

plants, the A6 group genes were the predominantly expressed 
HsfA genes. In particular, the A6e mRNA level in the roots 
in the 1.5 h heat treatment was ~25 times higher than that of 
any of the Al members. 

Subclass Bl and B2 members were also highly heat induc- 
ible with the 1.5 h treatment (Fig. 4). However, in contrast to 
B2 members, heat induction of the Bl group in the leaves was 
transient, and repeated heat stress for 3 d resulted in down- 
regulation of their transcript levels. The transcript levels 
of the B2 group during heat stress were the highest among 
TaHsf subclasses. The B4 group was not heat inducible and 
their transcript levels in the roots of plants with repeated heat 
stress were significantly lower than in the non-stress control 
plants. The CI group, which were expressed in the leaves 
under non-stress conditions at very high levels compared with 
other TaHsf groups, generally showed rapid down-regulation 
in the early heat treatment (1.5 h), particularly in the leaves, 
and their transcript levels gradually returned back to the con- 
trol levels with prolonged heat treatment. Most C2 members, 
which were not expressed in the leaves but were expressed at 
high levels in the control roots, were also down-regulated by 
heat stress. However, C2a and C2b expression was inducible 
in the leaves by >5-fold with a 1.5 h heat treatment, but was 
down-regulated after repeated heat stress for 3 d. 



Interestingly, heat induction of TaHsf genes was rapidly 
attenuated in the 5h treatment, and for many TaHsf genes 
further reduction was seen in samples repeatedly heat treated 
for 3 d (5h d~'). This heat response pattern was also observed 
in the TaHsp genes examined (Fig. 3B). 

A number of TaHsf subclasses are up-regulated during 
drought and salt stress 

Although the major role of Hsfs is known to be the regu- 
lation of heat-responsive genes involved in heat acclimati- 
zation, it was also of interest to see whether this family is 
involved in wheat responses to other major abiotic stresses 
such as drought and salt stress. Using the nucleotide 
sequences of these identified TaHsf members to search 
probe sets in the Affymetrix wheat genome array at the 
NetAffx Analysis Center of the Affymetrix website (http:// 
www.affymetrix.com/index.affx), 22 probe set IDs that 
have at least four matches out of 11 probes per probe set 
ID were identified (Table 1). However, eight probe set IDs 
cross-hybridize with more than one TaHsf gene, and two 
genes (C2a and C2g) are represented with each having two 
probe set IDs (Fig. 5). There are two publically deposited 
Affymetrix data sets, TA23 (http://www.plexdb.org; Aprile 
et al, 2009) and E-MEXP-971 (http:/www.ebi.ac/array- 
express/; Mott and Wang, 2007), for transcript profiling 
in wheat responses to drought and salt stresses. The TA23 
data set investigates expression changes in the flag leaves 
and glumes (two organs pooled) of hexaploid wheat (T. aes- 
tivum) and tetraploid wheat (T. turgidum subsp. durum) at 
the grain-filling stage during mild or severe drought stress 
in comparison with control plants. Among these Hsf probe 
sets, 13, 12, and 10 probe sets had detectable hybridization 
signals in the flag leaf and glume tissues of T. aestivum cv. 
Chinese Spring, the Chinese Spring-5A deletion line, and 
durum wheat, respectively. Of these detectable probe sets, the 
A2, A3, A6, CI, and C2 groups generally showed increased 
expression during drought stress. In contrast, the transcript 
levels of the A4 and Bl groups were down-regulated by 
drought. The drought-mediated expression level changes 
(<4-fold) of these Hsf genes were not as great as those seen 
in the short-term heat stress (Fig. 4), but were comparable 
with those during long-term heat stress (3 d with 5 h d~'). The 
transcript level changes of these TaHsf genes in the shoots 
and roots of T. aestivum plants after long-term salt stress 
from the E-MEXP-971 Affymetrix data set are shown in 
Fig. 6. In this data set, five genotypes were placed (each with 
three replicates) with the same treatment in the same group 
to increase the number of biological replicates for analysis. 
A total of 12 and eight probe sets with a significant increase 
in expression were observed in the salt-treated shoots and 
roots, respectively. The up-regulated TaHsf genes were A2, 
A6, A8, Bl, CI, and C2 groups in the shoot and A2, A4, 
A6, and CI groups in the root. Some of these increases in 
the expression levels in salt-stressed plants were remarkable. 
A > 10-fold increase in the A2d mRNA level was seen in the 
salt-stressed shoot, and a 10-fold increase in the A4 group 
expression levels was observed in the salt-stressed root. 
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Fig. 4. Expression changes of TaHsf members during heat stress. Values are means ±SD of three biological replicates and are 
expressed as apparent expression levels relative to a control gene TaRP15. Statistical significance in differences between control and 
heat-treated groups (36 °C for 1 .5h, 5h, or 3 d with 5h d~ 1 ) is indicated by a black (up-regulated) or blue (down-regulated) asterisk. 
(Note: the expression levels of some TaHsfs are too low to be seen in this comparative presentation, but statistical significance in the 
changes of their transcript levels in response to heat stress is shown.) 



TaHsfA2b, A4e, and C2a are transcriptional activators 
of heat shock protein genes 

To investigate the regulatory role of TaHsf transcription factors 
in controlling the expression of Hsp genes in wheat, transacti- 
vation assays were performed in wheat seedlings using the pro- 
moters of TaHspl 7 and TctHsp90.1-Al to drive the expression 
of gfp as reporter genes, and a constitutively expressed maize 
Ubil promoter to drive the expression of five TaHsf members 
from A2, A4, Bl, CI, and C2 subclasses (A2b, A4e, Bib, Clb, 
and C2a) as effector genes. TaHspl 7 and TaHsp90.1-Al pro- 
moter fragments were isolated based on the assembled sequences 
obtained through analysis of the wheat genome sequence data- 
base to extend the sequences of TaHspl 7 (CD909515 and 
CJ523462) and TaHsp90.1-Al (HX149366 and GQ240772) 
cDNA or ESTs. A control GUS reporter gene driven by the 
maize Ubil promoter was used for checking the transformation 
events across the samples. Quantitative real-time PCR analysis 



showed that the mRNA levels of TaHspl 7 and TaHsp90.1-Al 
were extremely low in wheat leaves and roots under non-stress 
conditions (Fig. 3). Green foci of the gfp reporter genes driven 
by the TaHspl7 or TaHsp90.1-Al promoter were rarely found 
in wheat seedling shoots and roots incubated at 22 °C, when 
they were bombarded into these seedlings without a TaHsf 
effector gene (Fig. 7). When the reporter gene-bombarded 
seedlings were incubated at 36 °C, many strong GFP foci were 
observed (Fig. 7), indicating that expression of the TaHspl 7 
or TaHsp90.1-Al promoter-driven reporter gene was induced 
by heat treatment. To confirm that the 22 °C samples bom- 
barded with the TaHspl 7 or TaHsp90.1-Al reporter gene were 
transformed, tissue sections where GFP foci could be found 
were selected for histochemical staining of GUS activity, and 
many strong GUS foci were observed in the samples incubated 
at 22 °C (Fig. 7), confirming that these gfp reporter genes were 
essentially not expressed in the shoots and roots of wheat 
seedlings at 22 °C. Introduction of these reporter genes with a 



548 I Xueefa/. 



5 
4.5 

4 
3.5 

3 
2.5 

2 

1.5 - 

1 
0.5 

0 



TT T 



CS Control 
I CS mild drought 
I CS Severe drought 

CS-5AL Control 
■ CS-5AL mild drought 
I CS-5AL Severe drought 

Creso Control 
I Creso mild drought 
I Creso Severe drought 



T T T T 



It 



m 



no* kj^ ^ q.^ k 1 ^ ^ ^ n 1 ^ <Y^ r- 1 ^ 

3" iv* *»* v v ,s> ,sS * # ,<v ,<w- <v c> ,c,V ,c,V .0" ,c,™ ,cy 



«jr or 



► _> *> *> 



cx.x' *y c.x' V' <^ N ' V Ci N ' *>V ^ c^' . V* .V* Ci^ \G t^' ."V' * 



Fig. 5. Drought responsiveness of TaHsf genes in the flag leaves and glumes of T. aestivum Chinese Spring (CS) and T. durum Creso. 
Raw data are derived from an Affymetrix wheat genome array data set at http://www.plexdb.org (accession # TA23; Aprile ef al. , 2009). 
Values are means ±SD of three biological replicates and relative expression levels within each genotype (stress groups versus control 
group and each control group was arbitrarily set as 1). 5AL, Chinese Spring 5A deletion line. Hybridization signal <20 is considered not 
detectable (nd). *P < 0.05. 



TaHsf effector gene (A2b, A4e, or C2a) resulted in induction of 
GFP foci in the shoots and roots of seedlings at 22 °C (Fig. 7; 
the activation pictures of the TaHsp90.1-Al reporter gene by 
HsfA4e or HsfC2a effector gene are not illustrated), indicating 
that TciHspl7 and TaHsp90.1-Al promoters were activated by 
a TaHsf transcriptional activator. GFP foci in the samples co- 
transformed with a TaHsf A2b effector gene were stronger than 
those with a TaHsf A4e effector gene. Interestingly, GFP foci 
in the samples co-transformed with a TaHsfCla effector gene 
were as strong as those with the TaHsf Alb effector gene. There 
was no induction of the TaHspl7- or TaHsp90.1-Al -driven 
reporter gene expression when TaHsfBlb or TaHsfClb was 
used as an effector gene (data not shown). 

TaHsfA2b binds to c'\s-acting elements in the 
promoters of heat shock proteins 

Transactivation of TaHsp 17 and TaHsp90,l-Al genes by the 
TaHsfA members suggests that these Hsp genes are their tar- 
get genes. To investigate further whether this transactivation 
is direct or indirect, whether the class A group was capable of 
directly binding to the ds-acting elements in the promoters 
of heat-inducible Hsp genes was studied using the TaHsfA2b 
protein. The available promoter sequences within 1.5 kb 
upstream of the translation start of four Hsp genes (TaHspl7, 
TaHsp90.1-Al, TaHsp26.6, and TaHsplOd) were analysed for 
the presence of HSEs. TaHsp26.6 and TaHsp70d promoter 
region sequences were also identified through analysis of the 
wheat genome sequence database to extend the sequences of 
TaHsp26.6 (AF097659) and TaHsp70d (AF005993) cDNAs 
(assembled TaHsp26.6 and TaHsp70d promoter region 
sequences are shown in Supplementary Fig. S5 at JXB online). 
A typical HSE ( 5 '-nG AAnnTTCnnG AAn or 5-nTCC- 
miGAAnnTTCn) was found only in the TaHsp90.1-Al 



promoter and was named TaHsp90.1El. Therefore, some 
sequences containing at least two nGAAn inverted repeats 
or three nGAAn inverted repeats with gaps or one nucleotide 
degeneracy in one repeat were retrieved for testing of their 
binding by TaHsfA2b protein in vitro using a CELD fusion 
system (Xue, 2002a). As shown in Fig. 8, TaHsf A2b binds 
to many of the sequences selected from these four promoters 
with varying affinity. At least one functional HSE was found 
in each of these four promoters. Some atypical HSE sequences 
(TaHsp26.6El and TaHsp70dE2) had even higher TaHsf A2b 
binding affinity than TaHsp90.1El. TaHsp90.1E3 with an 
additional three nucleotides in one spacer of three nGAAn/ 
nTTCn repeats had little detectable TaHsfA2b binding. As 
some functional TaHsp HSE sequences had one-nucleotide 
degeneracy, three synthetic sequences (GAAn2GTCn2GAA, 
TCCn2GCAn2TTC, and AGAAn2TTCT) were tested for 
TaHsf A2b binding activity. Either non-detectable or <10% 
of the TaHsp90.1El activity was observed in these synthetic 
oligonucleotides. These data indicate that single-nucleotide 
degeneracy in a perfect HSE sequence is allowed only in cer- 
tain positions and three nGAAn inverted repeats are required 
for effective binding by HsfA. In addition, comparison of 
TaHsp 17E1 and TaHsp 17E2 binding activity indicates that 
the sequence flanking the GAA or TTC core motif also affects 
Hsf binding. The binding analysis together with transactiva- 
tion assays implies that TaHspl7 and TaHsp90.1-Al are likely 
to be the direct target genes of TaHsfA2b. 

TaHsp90. 1E1 and TaHsp17E1 are functional heat 
shock elements and responsible for transactivation by 
TaHsf A2b and TaHsfC2a 

To investigate whether TaHsp90.1 El (GAAGCTTCGGGAA) 
and TaHsp 17E1 (GAACATTTTGGAA) are functional HSEs 
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Fig. 6. Expression changes of TaHsf genes in the shoot and 
roots of bread wheat plants in response to long-term salt stress. 
Raw data are derived from Affymetrix wheat genome array at 
http://www.ebi.ac.uk/arrayexpress/ (accession #: E-MEXP-971 ; 
Mott and Wang, 2007). Values are means ±SD of 1 5 biological 
replicates and are expressed as relative expression level (stress 
group versus control group and each control group was arbitrarily 
set as 1). The Affymetrix probe sets for these genes are in the 
same order as shown in Fig. 5. Hybridization signal <20 is 
considered not detectable (nd). *P < 0.01 . 



in vivo, TaHsp90.1-Al mutant promoters with truncation 
around the TaHsp90.1El region and TaHsp90.1-Al minimal 
promoters with an addition of TaHsp90.1El or TaHspl7El 
were tested (Fig. 9A). The gfp reporter gene driven by a 
short TaHsp90.1-Al promoter fragment (328 bp upstream 
of its translation start codon) (sHsp90gfp) was strongly acti- 
vated by heat treatment, while a deletion in TaHsp90.1El 
(AHSE90gfp) completely abolished its heat-inducible expres- 
sion (Fig. 9B). Addition of TaHsp90.1El or TaHspHEl to 
AHSE90gfp (the minimal TaHsp90.1-Al promoter) restored 
its heat-inducible expression (Fig. 9B). 



With these TaHsp90.1-Al mutant promoter reporter genes, 
it was investigated whether TaHsp90.1El and TaHspHEl are 
responsible for transactivation by TaHsfA2b and TaHsfC2a. 
As shown in Fig. 10, expression of the gfp reporter genes 
driven by a sHsp90, HSE90, or HSE17 promoter was strongly 
activated by both Ubil promoter-driven TaHsj A2b and 
TaHsfC2a effector genes without heat treatment. TaHsfA2b 
and TaHsfC2a were unable to activate the expression of a 
TaHsp90.1El deletion construct (AHSE90gfp) (Fig. 10), 
indicating that TaHsp90.1El and TaHspHEl are essential 
for transactivation by these two Hsf proteins. To confirm 
whether the transactivation of these reporter genes is specific 
to Hsf-HSE interaction, Ubil promoter-driven barley CBF1 
(UbilHvCBFl), which is a transcriptional activator of cold- 
inducible genes (Xue, 20026), was also tested. The HSE90gfp 
construct was not transactivated by the UbilHvCBFl effec- 
tor gene (Fig. 10). 

TaHsfC2a contains a functional AHA motif that is 
responsible for its activation activity 

Transactivation of TaHsp90. 1 E 1 - and TaHsp 1 7E 1 -containing 
promoters by TaHsfC2a is particularly interesting, as it is 
generally considered that HsfC class members do not contain 
a transcriptional activation domain. Therefore, it was inves- 
tigated whether an AHA-like sequence (LLLDGDFGNVSA 
FGPDAVDFAGFYTDDAFANAPVPVE) present in 
TaHsfC2a is responsible for its transactivation activity. Two 
TaHsfC2a mutant constructs (C2a-mCl and C2a-mC2) were 
made with mutation at the AHA-like region (Fig. 11 A). 
The UbilC2a-mC2 construct containing a replacement of 
three aromatic residues (FxxFY) with KxxSS and two acidic 
residues (DD) with KR abolished its transactivation activ- 
ity of the TciHspl7 reporter (Hspl7gfp) (Fig. 1 IB), while 
the mutation in the non-AHA residues of this AHA-like 
motif (UbilC2a-mCl) retained its transactivation activity. 
This experiment demonstrated that this AHA-like motif in 
TaHsfC2a is responsible for its transactivation activity. 



Discussion 

This study identified 56 Hsf genes from T. aestivum, which 
can be divided into three classes, A, B, and C. Although not 
all Hsf members have been identified in this study because 
the wheat genome sequence is still incomplete, analysis of the 
phylogentic tree of Hsf proteins from wheat and rice shows 
that the representatives of all rice subclasses have been identi- 
fied from wheat (Fig. 1). Phylogenetic analysis showed that 
TaHsf proteins are more closely clustered with rice Hsfs. 
Similar to rice, the wheat Hsf family does not contain A9 
and B3 subclasses as are present in Arabidopsis, but has C2 
subclass members. Phylogenetic analysis also indicated strong 
support for separate B and C classes of Hsfs, but class A Hsf 
protein members formed at least three smaller clusters. Two 
of the class A clusters appear to be more similar to class C 
Hsf proteins than to the third class A cluster. This observa- 
tion suggests that the A class is paraphyletic and that the C 
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Fig. 7. Transactivation of TaHsp promoter-driven reporter genes by TaHsfs. (A) Reporter and effector gene constructs. 
(B) Transactivation of the TaHspU or TaHsp90.1-A1 promoter-driven gfp reporter gene by a constitutively expressed TaHsfA2b, A4e, 
or C2a effector gene. The maize Ubil promoter-driven GUS reporter was co-bombarded with test constructs to check transformation 
events. GFP foci (green) indicate the expression of the TaHsp promoter-driven gfp reporter gene, and blue foci resulted from the 
expression of the co-introduced GUS reporter as indication that tissue sections were transformed with these constructs. The red 
background is the red fluorescence of shoot chlorophyll. The TaHsfA4e or TaHsfC2a effector gene also activated the expression of the 
Hsp90gfp reporter gene (data not shown). Co-bombardment of a TaHsfBlb or TaHsfClb effector gene with a Hsp17gfp or Hsp90gfp 
reporter gene did not produce GFP foci (data not shown). 



class is derived from a branch of the A class or that some 
subclasses of class A are really part of the C class. A similar 
observation was made by Scharf et al. (2012) when describ- 
ing the phylogenetic relationships across Hsfs from nine plant 
species including both monocots and eudicots. They also 
noted a paraphyletic grouping of the A class Hsf members, 
whereas the B and C classes of Hsfs were monophyletic. It 
should also be noted that the branching patterns seen in the 
current phylogenetic tree figures may be influenced by incom- 
plete data sets (missing wheat Hsfs). In comparison with 
Arabidopsis and tomato, monocot species (wheat and rice as 
well as Brachypodium distachyori) in general have fewer Al 
members, and more A2 and C members per genome (Table 1 ; 
Scharf et al, 2012). However, not much is known about the 
functional diversification of Hsfs in monocots to date. 

In this study, the detailed expression profiles of individual 
TaHsf members in six important organs and heat responses 



were provided. The subclass Al group members in wheat 
were found to be constitutively expressed at relatively high 
levels in all organs examined and were not up-regulated by 
heat stress, which is similar to HsfAl genes from other species 
such as Arabidopsis and rice (Hiibel and Schoffl, 1994; Mittal 
et al., 2009). However, there were many other constitutively 
expressed TaHsfA members (e.g. the A8 group, and some of 
the A2 and A6 group members) that were also expressed at 
relatively higher levels. A2 and A6 members were the main 
TaHsf members within class A that were up-regulated at 
a very high level during heat stress. These Hsf proteins are 
likely to be major transcription factors for up-regulation of 
Hsp genes during heat stress, as observed in wheat leaves and 
roots. Other class A genes that were up-regulated by heat 
were members of the A7 group and some A4 members, but 
they were expressed at a low level under heat stress, com- 
pared with the A2 and A6 groups. However, in terms of fold 
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Fig. 8. TaHsfA2b binding to elements present in the promoters of 
TaHsp17, TaHsp26.6, TaHp70d, and TaHsp90.1-A1 genes. 
(A) Relative binding activity (RBA) of TaHsfA2b to nGAAn or nTTCn 
sequences. Values are means ±SD of three assays and are relative 
to the binding activity of TaHsp90.1 E1 , which is arbitrarily set 
as 1 . GAA or TTC sequences are in bold. TaHsp sequences are 
selected from the available promoter region sequences within 
1 .5 kb upstream of the translation start of these TaHsp genes. 
Three artificial sequences (GAAn2GTCn2GAA, TCCn2GCAn2TCC, 
and AGAAn2TTCT) were also tested for tolerance of nucleotide 
sequences deviating from the perfect HSE (nGAAnnTCCnnGAA or 
nTCCnnGAAnnTCCn). (B) Illustration of TaHsfA2b binding activity 
as the 4-methylumbelliferyl group fluorescence produced through 
hydrolysis of methylumbelliferyl-p-D-cellobioside by TaHsfA2b- 
CELD that bound to immobilized HSE-containing oligonucleotides 
(TaHsp70dE1 and TaHsp90.1E1). Control is the oligonucleotide 
containing no HSE elements. 

induction, heat up-regulation of TaHsf A7 genes was very 
high. For example, the TaHsf A7 a mRNA level was increased 
by 127-fold in the leaves at the 1.5 h heat treatment (data not 
shown). In class B of the TaHsf family, the Bl and B2 groups 
were strongly up-regulated by heat; however, the expression 
levels of the Bl group gradually decreased after an initial 
strong induction. In class C, the CI group were constitutively 
expressed at a high level in most organs. In young leaves, 
the expression levels of the CI group were even higher than 
those of the Al group The CI group members were generally 
down-regulated by heat stress, and their heat response pattern 
appears to be in the opposite direction to that of the heat-up- 
regulated A2 and A6 group genes. The C2 group members 
were predominantly expressed in the embryos and roots, and 
were also generally down-regulated by heat stress, except C2a 
and C2b that were up-regulated in the leaves by a short-term 
heat treatment. Interestingly, marked up-regulation of HsfC2 
genes from rice by oxidative and heat stresses or a combina- 
tion of these stresses has been reported (Mittal et al, 2009, 
2012). 

From analysis of TaHsf expression profiles, two interest- 
ing features are worth mentioning here. First, a number of 
Hsf members other than the Al group were constitutively 
expressed at a level similar to or higher than the Al group, 
presuming that they are generally present in an inactive form 
as is the case for HsfAl proteins known in other species 
(Hahn et al, 2011; Scharf et al, 2012). Secondly, in contrast 
to what has been observed in sunflower, Arabidopsis, and 
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Fig. 9. The heat-induced expression of gfp reporter genes driven 
by a truncated TaHsp90. 1-A1 promoter or a minimal promoter 
with an addition of TaHsp90.1E1 (GAAGCTTCGGGAA) or 
TaHsp17E1 (GAACATTTTGGAA). (A) Reporter gene constructs. 
A short TaHsp90. 1 -A 1 promoter (sHsp90) has a 328 bp fragment 
upstream of the translational start codon, which contains the 
TaHsp90.1 E1 element. The sequence upstream the TATA box 
of the TaHsp90. 1-A1 promoter is shown. The AHsp90 promoter 
starts immediately downstream of the TaHsp90.1 E1 element. 
HSE90 and HSE1 7 constructs contain three TaHsp90.1 E1 or 
TaHspl 7E1 repeats, which are added immediately upstream of 
the AHsp90 promoter (a TaHsp90.1-Al minimal promoter). 
(B) Transient expression assays of reporter genes. The UbilGUS 
reporter gene was also co-introduced, and GUS foci are illustrated 
when gfp expression was essentially undetectable. 

rice (Almoguera et al., 2002; Kotak et al., 2007; Chauhan 
et al., 201 la), no seed-specific Hsf members were identified in 
wheat. However, three A6 members (A6c, A6d, and A6e) and 
A5b were predominantly expressed in endosperm under non- 
stress conditions. In addition, four A2 members (A2b, A2c, 
A2e, and A2f) were also expressed in the endosperm of non- 
stressed plants at a relatively high level. This suggests that 
the endosperm is rich in TaHsfA proteins under non-stress 
conditions, with the assumption that there is no major dif- 
ference at the translational level among TaHsfA subclasses. 
The endosperm is an important organ for wheat productiv- 
ity. Some enzymes involved in starch synthesis in several crop 
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Fig. 10. TaHsfA2b and TaHsfC2a transactivate the gfp reporter 
driven by the TaHsp90. 1-A1 minimal promoter containing 
TaHsp90.1E1 (HSE90) orTaHsp17E1 (HSE1 7). Constructs 
used are shown in Fig. 9. The UbilGUS reporter gene was also 
co-introduced, and GUS foci are illustrated when gfp expression 
was essentially undetectable. Ubil promoter-driven HvCBFI , 
which is a transcriptional activator for cold-inducible genes, was 
used as a negative control. 

species are known to be thermolabile, such as soluble starch 
synthase from wheat endosperm (Keeling et al. , 1993) and pea 
embryos (Denyer and Smith, 1992) and ADP-glucose pyroph- 
osphorylase from maize endosperm (Greene and Hannah, 
1998). Biochemical analysis has shown that the crude extract 
of soluble starch synthase from wheat endosperm loses 
its enzyme activity significantly with a 15min treatment at 
30 °C, and a decrease in the enzyme activity is noticeable even 
at 25 °C with a 2h treatment (Keeling et al, 1993). Wheat 
crops in Australia and other warm climate countries often 
encounter high temperatures above 30 °C at the grain-filling 
stage, which results in a significant reduction in grain yield 
and quality (Skylas et al, 2002). The A6 and A2 Hsf genes 
that were highly expressed in the endosperm are likely to be 
important Hsf members for keeping the expression of Hsp 
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Fig. 1 1 . Loss of the transactivation activity of TaHsfC2a by 
mutation of hydrophobic and acidic residues in its C-terminal 
domain. (A) TaHsfC2a and its mutant constructs. The substituted 
amino residues in the mutants (C2a-mC1 and C2a-mC2) are 
in red. (B) Transactivation of the TaHspl 7 promoter-driven gfp 
reporter (Hsp1 7gfp) by TaHsfC2a and its mutants. The Ubil GUS 
reporter gene was also co-introduced, and GUS foci are illustrated 
when gfp expression was essentially undetectable. 

genes at a high level. Indeed, four out of the five TaHsp genes 
examined were also expressed at a high level in the endosperm 
under non-stress conditions, while their mRNA levels in the 
vegetative organs were either not detectable or at a much 
lower level. It is tempting to speculate that this expression 
pattern of heat-inducible TaHsp genes may have an advan- 
tage in protecting this important organ against heat stress. 
Many Hsps are known to function as molecular chaperones 
and are capable of stabilizing thermolabile proteins against 
heat denaturation (Grover et al, 2013; Waters, 2013). For 
example, a mitochondrial small Hsp can improve the electron 
transport activity of the mitochondrial NADH: ubiquinone 
oxidoreductase complex I in plants at high temperatures 
(Downs and Heckathorn, 1998). Similarly, a small Hsp from 
the chloroplast protects the chloroplast photosystem II elec- 
tron transport activity from heat inactivation (Hechathorn 
et al., 1998). 

Besides their role in wheat adaptation to heat stress, 
a number of TaHsf groups (A2, A3, A4, A6, Bl, and C2) 
were up-regulated in the flag leaves and glumes of wheat at 
the grain-filling stage by drought stress. The up-regulation of 



these genes as well as the A8 and CI groups was also seen in 
the shoots of salt-stressed plants. A similar salt up-regulation 
of these genes was observed in the roots, with the exception 
of the A8 and C2 groups. These data indicate that some of 
TaHsf members may have a regulatory role in wheat adapta- 
tion to drought and salt stress. Intriguingly, an examination 
of the expression levels of TaHspl6.9, TaHspl7, TaHsp26.6, 
and TaHsp90.1 genes in drought- and salt-stressed wheat 
plants in the Affymetrix wheat genome array data sets (TA23 
and E-MEXP-971) did not reveal up-regulation of these 
TaHsp genes by drought and salt stresses (data not shown). 
This suggests that drought- or salt-induced up-regulation of 
these TaHsfs has no impact on the expression of these TaHsp 
genes. Possible explanations for this observation include that 
TaHsf proteins are inactive, or that Hsp genes are in a strongly 
repressed state in these salt- and drought-stressed organs. 
Whether these TaHsfs are able to regulate the expression of 
other stress protection genes during salt and drought stresses 
requires future investigation. However, it has been shown that 
overexpression of some HsfA genes, such as AtHsfA2 and 
HaHsfA9, can result not only in thermotolerance but also 
in salt and dehydration stress tolerance (Ogawa et al. 2007; 
Prieto-Dapena et al., 2008), implying that some drought- and 
salt-inducible TaHsf members may have a similar role. 

In this study, some functional HSEs present in the promot- 
ers of TaHspl7, TaHsp26.6, TaHsp70d, and TaHsp90.1-Al 
genes, which can be strongly bound by TaHsfA2b, were identi- 
fied. Most of these functional HSEs in the promoters of these 
wheat Hsp genes contain imperfect HSE sequences, some of 
which had TaHsfA2b binding affinity even higher than that of 
the perfect HSE in the TaHsp90.1-Al promoter. It was shown 
that TaHsfA2b and TaHsfA4e were transcriptional activators 
and capable of activating the expression of TaHspl7- and 
TaHsp90.1-Al promoter-driven reporter genes in the leaves 
and roots. It appears that TaHsfA2b is a more potent transcrip- 
tional activator in transactivation of the TaHsp 1 7 promoter- 
driven reporter gene than TaHsf A4e. However, transactivation 
analysis showed that TaHsfBlb and TaHsfClb had no ability 
to activate TaHsp 17 and TaHsp90.1-Al expression. As no 
direct evidence in previous studies has shown that Bl and CI 
members are transcriptional activators of Hsp genes, it is likely 
that TaHsfBl and TaHsfCl proteins are not responsible for 
up-regulation of these Hsp genes in wheat, although TaHsfBl 
genes were up-regulated 5- to 10-fold in the leaves and roots at 
the first 1.5 h heat treatment. In heat-stressed tomato, HsfBl 
acts as a transcription co-activator, cooperating with class 
A Hsfs (Bharti et al., 2004), while Arabidopsis HsfBl func- 
tions as a repressor (Czarnecka-Verner et al., 2004). However, 
rice HsfCla and HsfClb can potentially serve as transcrip- 
tional activators, as the transactivation activity of these two 
HsfCl proteins has been observed in yeast cells by fusion with 
the Gal4 DNA-binding domain (Mittal et al. , 201 1). 

Promoter truncation and mutation analyses showed 
that the TaHsp90.1El (GAAGCTTCGGGAA) of the 
TaHsp90.1-Al promoter is the HSE responsible for its 
heat-induced expression by Hsf proteins. An atypical HSE 
identified in the TaHspl7 promoter, TaHsp 17E1 (GAACA 
TTTTGGAA), is also a functional HSE in vivo. The addition 
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of either TaHsp90.1El or TaHsp 17E1 to a minimal promoter 
derived from the TaHsp90.1-Al promoter resulted in the 
heat-inducible expression and transactivation by TaHsfA2b 
or TaHsfC2a at 22 °C. These data implicate that the presence 
of either TaHsp90.1El or TaHspl7El is sufficient for heat- 
mediated expression through interaction with Hsf proteins. 

Most interestingly, it was observed that TaHsfC2a acted 
as a transcriptional activator and was capable of transac- 
tivating the expression of the TaHspl7 or TaHsp90.1-Al 
promoter-driven reporter gene. No previous studies to 
date have shown that a class C Hsf protein can transacti- 
vate the Hsp genes. The C-terminal end of TaHsfC2a con- 
tains a motif, LLLDGDFGNVSAFGPDAVDFAGFYTDD 
AFANAPVPVE, which resembles an AHA motif (Nover 
and Scharf, 1997; Kotak et al., 2004). The mutation of some 
aromatic and acidic residues in this motif resulted in the abol- 
ishment of TaHsfC2a transactivation activity, suggesting that 
this AHA-like motif is responsible for its transcriptional acti- 
vator property. TaHsfC2a was predominantly expressed in 
the seeds, particularly in the embryo, but was also heat induc- 
ible in the leaves with a short-term heat treatment (a 5-fold 
increase). Up-regulation of rice HsfCla by a short-term heat 
treatment (10-30 min) has also been shown in rice seedlings 
(Mittal et al, 2009). An attempt was also made to investigate 
whether TaHsfC2a is capable of binding to the functional 
HSEs present in the TaHsp promoters in vitro, but the experi- 
ment failed technically due to a very low level of TaHsfC2a- 
CelD fusion protein produced in Escherichia coli. However, 
the transactivation of TaHspl7 and TaHsp90.1-Al genes by 
TaHsfC2a appears to be through its specific interaction with 
TaHspl7El and TaHsp90.1El elements, as the promoter 
mutagenesis experiment demonstrated that the transactiva- 
tion of gfp reporter genes by TaHsfC2a relied on the presence 
of either TaHsp90.1El or TaHsp 17E1. Based on its expres- 
sion pattern and ability to activate Hsp genes, TaHsfC2a 
probably plays a role in regulation of Hsp expression in the 
embryos and during early heat stress. 

Another intriguing point raised by this study is that 
transactivation of Hsp genes by TaHsfA2b, TaHsfA4e, and 
TaHsfC2a does not require heat treatment, indicating that 
TaHsf proteins are functional at 22 °C. Presumably, they 
are able to form functional trimers under non-stress condi- 
tions. The question is why heat-inducible Hsp genes such as 
TaHspU are not expressed in the leaves and roots of wheat 
under normal conditions, as some A2 (A2a, A2d, and A2f) 
and C2 (C2d and C2g) members were expressed in these 
organs at a level similar to or higher than that of TaHsf A 1 
members. It has been shown that some repressors (e.g. Hsp70 
and Hsp90) bind to Hsf protein to keep them in an inactive 
monomer form (Kim and Schoffl, 2002; Yamada et al. 2007; 
Hahn et al, 2011). At present, a model of Hsf activation 
derived from tomato plants is that HsfAl is a master regulator 
of the heat response (Mishra et al. , 2002), and is in an inac- 
tive state through complexing with Hsp90 and Hsp70 (repres- 
sors) in non-stress conditions (Hahn et al. , 201 1 ; Scharf et al. , 
2012). In Arabidopsis, there is no single master Hsf, and three 
Al members together serve as heat response-triggering factors 
(Liu et al, 2011; Liu and Charng, 2012). Upon heat stress, 
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Hsp90 and Hsp70 are released from HsfAl to form a func- 
tional trimer, which subsequently activates expression of other 
HsfA members. However, constitutively expressed HsfA genes 
in wheat include some A2, A6, and A8 members in addition 
to Al members. This implies that some repressors can keep 
all these constitutively expressed TaHsf transcriptional activa- 
tors in an inactive form in wheat. An alternative hypothesis 
is that some repressors (e.g. transcriptional repressors) other 
than Hsp90 and Hsp70 may interfere with Hsf transcriptional 
activity at the Hsf protein and Hsp promoter assembly point 
in wheat to keep Hsp genes in a repressed state under non- 
stress conditions. These hypotheses suggest that Hsp gene 
activation by the Hsf proteins in wheat plant cells is likely to 
be dependent on the relative amounts of TaHsf proteins and 
their repressors. When TaHsf protein amounts are in a sub- 
stantial excess compared with the amounts of their repressors 
such as Hsf overexpression in the transactivation assays or 
heat stress, Hsp genes become activated. This may also apply 
to the high level expression of some heat-inducible Hsp genes 
(e.g. TaHspl7, TaHsp26.6, TaHsp70d, and TciHsp90.1Al) in 
wheat endosperms and embryos under non-stress conditions. 

It will be interesting to investigate whether genetic vari- 
ability in Hsp levels and thermotolerance in wheat detected 
in early studies (Zivy, 1987; Krishnan et al, 1989; Vierling 
and Nguyen, 1992; Skylas et al, 2002) is contributed by dif- 
ferences in expression of some critical regulators — TaHsfs. 
Possibly, some TaHsf genes could be important quantita- 
tive trait loci for determining thermotolerance. A substantial 
amount of research is still required to unlock the mystery of a 
Hsf-Hsp regulation model and the molecular mechanisms of 
heat adaptation in this important crop species — wheat. This 
will be an important research issue to address in the future, 
to assist in ensuring an adequate food supply in the predicted 
global warming scenario. 
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